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We present a study of the decay B~ —^D^Qp-^K~ and its charge conjugate, where is re¬ 

constructed in CP-even, CP-odd, and non-CP flavor eigenstates, based on a sample of 232 million 
T(4S)—decays collected with the BABAR detector at the PEP-II e''"e“ storage ring. We measure 
the partial-rate charge asymmetries Acp± and the ratios Rcp± of the B^D^K decay branching 
fractions as measured in CP± and non-CP decays: Acp+ = 0.35 ± 0.13(stat) ± 0.04(syst), 
Acp- = —0.06 ± 0.13(stat) ± 0.04(syst), Rcp+ ~ 0.90 db 0.12(stat) db 0.04(syst), Rcp- = 0.86 ± 
O.lO(stat) dz 0.05(syst). 

PACS numbers: 11.30.Er,13.25.Hw,14.40.Nd 


A stringent test of the flavor and CP sector of the Stan¬ 
dard Model can be obtained from the measurements, in B 
meson decays, of the sides and angles of the unitarity tri¬ 
angle, which are related to the elements of the Cabibbo- 
Kobayashi-Maskawa matrix V. A theoretically clean 
measurement of the angle 7 = a,rg{—VudV*i,/VcdV*f,) can 
be obtained from the study of B ~decays Q 
by exploiting the interference between the b cus and 
b ucs decay amplitudes H i Among the pro¬ 
posed methods, the one originally suggested by Gronau, 
London, and Wyler (GLW) exploits the interference be¬ 
tween B~^D^K~ and B~K~ when the and 
mesons decay to the same CP eigenstate. 

The results of the GLW analyses are usually expressed 
in terms of the ratios Rcp± of charge-averaged partial 
rates and of the partial-rate charge asymmetries Acp±, 


Rcp± = 


Acp± = 


T{B-^D°Cp^K-) + TiB+^D°cp^K+) 
[r{B-^D0K-) + r{B+-^D°K+)] /2 ’ 
r{B-^D°CP^K-) - r{B+^D°cp^K+) 

TiB-^D°cp^K-) + riB+^D°CP±K+) ’ ^ ’ 


Here, ^ D^)lV2 are the CP eigenstates 

of the neutral D meson system, and we have followed 
the notation used in 0. Neglecting D^—D^ mixing 
the observables Rcp± and Acp± are related to the an¬ 
gle 7 , the magnitude r of the ratio of the amplitudes 
for the processes B~K~ and B~^D^K~, and the 
relative strong phase S between these two amplitudes, 
through the relations Rcp± = 1 -I- ± 2r cos S cos 7 and 

Acp± = ±2r sin(5sin7/i?cp± Q- Theoretical expecta¬ 
tions for r are in the range « 0.1 — 0.2 Hi , in agree¬ 
ment with the 90% G.L. upper limits on r set by BABAR 
(r < 0.23) and Belle (r < 0.1^ through the study of 
B~^DK~,D—fK^TT~ decays [3. 

In this paper we present the measurements of Rcp± 
and Acp±. The ratios Rcp± are computed using the 
relations Rcp± — R±/R, where the quantities R and R± 


are defined as: 


B{B-^D°K-) + B{B+^D°K+) 
BiB-^DOn-) + B{B+^D0Tr+) ’ 
B{B-^D°,p±K~) + B{B+^D°cp±K+) 
BiB-^D°CP^n-) + BiB+^D°cp±^+) 


(3) 

(4) 


Several systematic uncertainties cancel out in the mea¬ 
surement of these double ratios. We also express the 
CP-sensitive observables in terms of three independent 
quantities: 


x± 


r 


2 


_ Pcp+(1 T Acp+) — Pcp-(1 T Acp-) 

~ 4 

2 , 2 Rcp+ + Rcp- — 2 
— x± + y± — 7, ! 


(5) 

( 6 ) 


where x± = r cos((5 ± 7 ) and y± = r sin((5 ± 7 ) are the 
same CP parameters as were measured by the BABAR 
Collaboration with B~^DK~, D^KgTT~TT~^ decays 0. 
This choice allows the results of the two measurements 
to be expressed in a consistent manner. 

The measurements use a sample of 232 million T(4S') 
decays into BB pairs, corresponding to an integrated lu¬ 
minosity of 211 fb~^, collected with the BABAR detec¬ 
tor at the PEP-II asymmetric-energy B factory. Since 
the BABAR detector is described in detail elsewhere H, 
only the components that are crucial to this analysis are 
summarized here. Charged-particle tracking is provided 
by a five-layer silicon vertex tracker (SVT) and a 40- 
layer drift chamber (DCH). For charged-particle identi¬ 
fication, ionization energy loss in the DCH and SVT, 
and Cherenkov radiation detected in a ring-imaging de¬ 
vice (DIRC) are used. Photons are identified by the 
electromagnetic calorimeter (EMC), consisting of 6580 
thallium-doped Csl crystals. These systems are mounted 
inside a 1.5-T solenoidal superconducting magnet. We 
use the GEANT [l^ software to simulate interactions of 
particles traversing the detector, taking into account the 
varying accelerator and detector conditions. 

We reconstruct B~^D^h~ decays, where the prompt 
track h~ is a kaon or a pion. candidates are recon¬ 
structed in the CP-even eigenstates 7r“7r+ and K~K'^ 
{D^P_^_), in the CP-odd eigenstates Kgcj) and KgUj 
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{Dqp_), and in the non-CP, flavor eigenstate 
(j) candidates are reconstructed in the K~K'^ channel 
and u candidates in the 7r“7r+7r'^ channel. We optimize 
our event selection to minimize the statistical error on 
the B~^D^qp~^K~ signal yield, determined for each 
decay channel using simulated signal and background 
events. 

The prompt particle h is required to have a momen¬ 
tum greater than 1.4 GeV/c and the number of pho¬ 
tons associated to its Cherenkov ring is required to be 
greater than four to improve the quality of the recon¬ 
struction. We reject a candidate track if its Cherenkov 
angle does not agree within four standard deviations (cr) 
with either the pion or kaon hypothesis, or if it is iden¬ 
tified as an electron by the DCH and the EMC. Particle 
identihcation (PID) information from the drift chamber 
and, when available, from the DIRC, must be consistent 
with the kaon hypothesis for the K meson candidate in 
D^^K~Tr^, , and (j)^K~K^ decays and 

with the pion hypothesis for the meson candidates in 
P)°^7r“7r+ and decays. 

Neutral pions are reconstructed by combining pairs 
of photon candidates with energy deposits larger than 
70 MeV that are not matched to charged tracks. The 
77 invariant mass is required to be in the range 115-150 
MeV/c^ and the total 7r° energy must be greater than 
200 MeV. When 7r°’s are combined with other particles 
to form composite particles, the mass is constrained to 
the nominal mass [llj. 

Neutral kaons are reconstructed from pairs of op¬ 
positely charged tracks with invariant mass within 
7.8 MeV/c^ 3cr) of the nominal mass. We also re¬ 
quire that the ratio between the flight length in the plane 
transverse to the beam direction and its error be greater 
than 2. The cj) mesons are reconstructed from two op¬ 
positely charged kaons with invariant mass in the range 
1.008 < M{K^K~) < 1.032 GeV/c^. We also require 
|cos0hei(0)| > 0.4, where 0hei(0) is the angle between 
the flight direction of one of the 4> daughters and the 
flight direction, in the (j) rest frame. The uj mesons are 
reconstructed from Tr+Tr^Tr*^ combinations with invariant 
mass in the range 0.763 < M(7r+7r“7r°) < 0.799 GeV/c^. 
We define 0^ as the angle between the normal to the ui 
decay plane and the momentum in the co rest frame, 
and 07 r 7 r as the angle between the flight direction of one of 
the three pions in the u rest frame and the flight direction 
of one of the other two pions in their center-of-mass (CM) 
frame. The quantities cos9n and cos0,n.7r follow cos^0jv 
and sin^ 9^^^^ distributions for the signal and are almost 
flat for wrongly reconstructed or false uj candidates. We 
require the product cos^ 9n sin^ > 0.08. The invari¬ 
ant mass of a candidate, M{D^), must be within 2.5tT 
of the mean fitted mass, with resolution a ranging from 
4 to 20MeV/c^ depending on the decay mode. For 
Zl°^7r“7r+, the invariant mass of the (/i“7r+) system, 
where 7r+ is the pion from and h~ is the prompt track 


from B taken with the kaon mass hypothesis, must 
be greater than 1.9 GeV/c^ to reject background from 
B-^D°Tr-,D°^K-Tr+ and B-^K*°Tr-, K*°^K-Tr+ 
decays. When reconstructing B mesons, for all decay 
channels the candidate invariant mass is constrained 
to the nominal mass im. 

We reconstruct B meson candidates by combining a 
11° candidate with a track h. For the 
mode, the charge of the track h must match that of 
the kaon from the 11° meson decay. We select B me¬ 
son candidates using the beam-energy-substituted mass 
Toes = V (^0^/2 + Po • Pb) V^o “ p'b and the energy 
difference AE = Ep — E^j^^ where the subscripts 0 and 
B refer to the initial e+e” system and the B candidate 
respectively, and the asterisk denotes the CM (T(4S')) 
frame. The toes distributions for B~^D^h~ signals are 
Gaussian functions centered at the B mass with a reso¬ 
lution of 2.6MeV/c^, which do not depend on the decay 
mode or on the nature of the prompt track. In contrast, 
the AE distributions depend on the mass assigned to 
the prompt track and on the Zl° momentum resolution. 
We evaluate AE with the kaon mass hypothesis so that 
the distributions are Gaussian and centered near zero for 
B~^D^K~ events and shifted by approximately 50 MeV 
for B~^D^tt~ events. The B~^D^K~ AE resolution 
is about 17 MeV for all the 11° decay modes. All B can¬ 
didates are selected with toes within 3a of the peak value 
and with AE in the range —0.16 < AE < 0.23 GeV. 

To reduce background from continuum production 
of light quarks, we construct a linear Fisher discrimi¬ 
nant na based on the following quantities: (i) Lq = 
Pi and L 2 = J2i Pi 9i, evaluated in the CM frame, 
where pi is the momentum, and 9i is the angle with re¬ 
spect to the thrust axis of the B candidate of charged 
tracks and neutral clusters not used to reconstruct the 
B] (ii) I cos0t|, where 9t is the angle between the thrust 
axes of the B candidate and of the remaining tracks and 
clusters, evaluated in the CM frame; (iii) | cos 0 b |, where 
9b is the polar angle of the B candidate in the CM frame. 

For events with multiple B~^h~candidates (1-7% 
of the selected events, depending on the 11° decay mode), 
we choose that with the smallest formed from the dif¬ 
ferences of the measured and true masses of the candidate 
B, Zl°, 7r° (only for D°^K^Tr°, 0 

and to (D^^KgUj), compared to the appropriate recon¬ 
structed mass resolutions. The total reconstruction ef¬ 
ficiencies, based on simulated signal events, are 39% 
{K-TT+), 31% {K-K+), 30% (7r-7r+), 17% 20% 

and 7% {K°u;). 

The main contributions to the background from BB 
events come from the processes B^D*h {h = Tr,K), 
B~^D^p~, mis-reconstructed B~^D'^h~, and from 
charmless B decays to the same final state as the sig¬ 
nal: for instance, the process B~ —fK~K'^K~ is a back¬ 
ground for B~^D^K~, . These charmless 

backgrounds have similar AE and toes distribution as 
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the D^K~ signal and we call them “peaking BB back¬ 
grounds” . 

For each decay mode an extended unbinned max¬ 
imum likelihood fit to the selected data events deter¬ 
mines yields for two signal channels, B~ —>D^Tr~ and 
B~^D^K~, and four kinds of backgrounds: candidates 
selected either from continuum or from BB events, in 
which the prompt track is either a pion or a kaon. The 
fit uses as input AE and a particle identification proba¬ 
bility for the prompt track based on the Cherenkov angle 
9c, the momentum p, and the polar angle 6 of the track. 

The extended likelihood function £ is defined as 


N 


£ = exp - ^ Tii 


2=1 


J = 1 




. 2=1 


(7) 


where N is the total number of observed events and rii 
is the yield of the event category. The six functions 
Vi{xj-,di) are the probability density functions (PDFs) 
for the variables Xj, given the set of parameters di. They 
are evaluated as a product Vi = Vii{AE) x V 2 i{dc)- 

The AE distribution for B~K~ signal events is 
parametrized with a Gaussian function. The AE dis¬ 
tribution for >£)°7r“ is parametrized with the same 
Gaussian function used for B~ K~ with an addi¬ 
tional shift, computed event by event as a function of the 
prompt track momentum, arising from the wrong mass 
assignment to the prompt track. The offset and width 
of the Gaussian functions are determined from data to¬ 
gether with the yields. 

The AE distribution for the continuum background 
is parametrized with a linear function whose slope is 
determined from off-resonance data. The AE distribu¬ 
tion for the non-peaking BB background is empirically 
parametrized with the sum of a Gaussian function and 
an exponential function when the prompt track is a pion, 
and with an exponential function when the prompt track 
is a kaon. The parameters are determined from simulated 
events. The AE distribution for the peaking charmless 
BB background is parametrized with the same Gaussian 
function used for the B~ K~signal. The yield of the 
BB peaking background is estimated from the sidebands 
of the invariant mass distribution and fixed in the fit. 

The parametrization of the particle identification PDF 
is performed by fitting with two Gaussian functions the 
background-subtracted distribution of the difference be¬ 
tween the reconstructed and expected Cherenkov angles 
of kaon and pion samples. The parametrization is per¬ 
formed as a function of the momentum and polar angle 
of the track. Pions and kaons are selected from a pure 
£)*+ ^ £>°7r+, >i£“7r+ control sample. 

The results of the fit are summarized in Table n Fig¬ 
ure n shows the distributions of AE for the CP-f 

and CP— modes after enhancing the B^D^K purity 
by requiring that the prompt track be consistent with 
the kaon hypothesis. The total PDF, normalized by the 


fitted signal and background yields, integrated over the 
Cherenkov angle variable and modified to take into ac¬ 
count the tighter selection criteria, is overlaid in the fig¬ 
ure. 


TABLE I: Yields from the maximum likelihood fit. 


D'^ mode 

N{D'k+) 

N{DTr-) N{DK+) N{DK-) 

A'-7r+ 

8151 ±95 

7899 ± 93 

649 ± 29 

611 ± 28 

K-K+ 

705 ± 28 

690 ± 28 

26 ±9 

70± 10 

7r“7r’’" 

256 ± 18 

219 ± 17 

18 ± 7 

17 ±7 


707 ± 29 

677 ± 29 

39 ±9 

42 ± 9 

kU 

176 ± 14 

157 ± 13 

15 ±5 

13 ±4 

K°sio 

235 ± 17 

230 ± 17 

25 ±7 

14 ±6 



FIG. 1: Distributions of AE for events enhanced in 
B-^D°K- signal. Top: R-R-, D°^K-tt+- middle: 
B~^D%p^K~ ■, bottom: B~^Dqp_K~ . Solid curves rep¬ 
resent projections of the maximum likelihood fit; dashed, 
dashed-dotted and dotted curves represent the B^D^K, 
B^D^tv and background contributions. 


The ratios Rcp± are computed for the five CP modes 
using the relations in Eqs. 0 and ©■ A number of 
systematic uncertainties, such as the uncertainty from 
the tracking efficiency and the uncertainty on the de¬ 
cay branching fractions, cancel out in the measurement 
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of the double ratio. The relations Rcp± = R±/R hold 
under the assumption that the magnitude Ttt of the ra¬ 
tio of the amplitudes of the B~^D^Tr~ and B~ 
processes can be neglected @ (r^ ~ r < 0 . 012 , 
where A « 0.22 [m is the sine of the Cabibbo angle). 
This assumption is considered further when we discuss 
the systematic uncertainties. The quantities R±/R are 
computed from the ratios of the B^DK and B^Dtt 
yields in Table HJ scaled by correction factors taking into 
account small differences in the selection efficiency be¬ 
tween B^DK and B^Dtt. These correction factors 
are evaluated from simulated events and range between 
0.982 ± 0.018 and 1.020 ± 0.031 depending on the 
decay mode. The results for the CP-even and CP-odd 
combinations are listed in Table HD 

The partial-rate charge asymmetries Acp± are calcu¬ 
lated from the measured yields of positive and negative 
B^DK decays in Tabled The results for the CP-even 
and CP-odd combinations are reported in Table m 

TABLE II: Measured ratios Rcp± and Acp± for CP-even 
and CP-odd D decay modes. The first error is statistical, the 
second is systematic. Rcp- and Acp- are corrected for the 
CP-even dilution described in the text. 


mode 

Rcp 

Acp 

CP-t 

0.90 ±0.12 ±0.04 

0.35 ±0.13 ±0.04 

CP- 

0.86 ±0.10 ±0.05 

-0.06 ±0.13 ±0.04 


In the case of D^^K^cj), K~, and 

a>^7r'''7r“7r°, the values of Rcp- and Acp- quoted in Ta¬ 
ble HD are obtained after correcting the measured values 
to take into account the dilution from a CP-even back¬ 
ground arising from B~^D^h~, D^^Kg{K~K~^)non-!i> 
and D°^Kg{TT~Tr~^TT^)non-uj decays. For the Kgp chan¬ 
nel we exploit the investigation performed by BABAR 
of the K~ Dalitz plot to estimate the 

level of the CP-even background (0.160 ± 0.006 relative 
to the signal) and the corresponding Rcp- and 

Acp- dilution. For the KgW channel there is little in¬ 
formation on this background. We estimate the amount 
of P°^Pi 5 ( 7 r+ 7 r“ 7 r°)non-w background (0.25 ± 0.05 rel¬ 
ative to the K'^uj signal) from the cos9n distribution 
of B~^D^tt~, candidates, and assume 

the CP-even content of this background to be (50±29)%. 

Systematic uncertainties in the ratios Rcp± and in the 
CP asymmetries Acp± are listed in Table EH They 
arise both from the uncertainties on the signal yields, 
extracted through the maximum likelihood fit, and from 
the assumptions used to compute Rcp± and Acp±. The 
correlations between the different sources of systematic 
errors, when non-negligible, are considered when combin¬ 
ing the two CP-even or the three CP-odd modes. 

The uncertainties on the fitted signal yields are due 
to the imperfect knowledge of the AP and PID PDFs 


TABLE III: Systematic uncertainties on the observables 
Rcp± and Acp± after combination of the two CP-even and 
the three CP-odd decay modes. 


Source 

ARcp+ 

(%) 

ARcp- 

(%) 

AAcp+ 

(%) 

AAcp- 

(%) 

bkg. AP PDF 

1.3 

1.1 

1.1 

0.4 

PID PDF 

0.1 

0.1 

0.2 

0.2 

peaking bkg. yields 

3.0 

4.2 

2.6 

2.2 

opposite-CP bkg. 

- 

1.3 

- 

1.0 

detector charge asym. 

1.0 

1.1 

2.7 

2.7 

r-K 

2.2 

2.1 

- 

- 

Total 

4.1 

5.1 

3.9 

3.7 


and of the peaking background yields, and are evaluated 
by varying the parameters of the PDFs and the peak¬ 
ing background yields by ±lcr and taking the difference 
in the signal yields. The uncertainties in the branching 
fractions used in the simulation of the B decays that 
contribute to the BB background are also taken into 
account. The yields of the BB and continuum back¬ 
grounds found in data are consistent with what is ex¬ 
pected from the simulation. In the Kgp and KgUj chan¬ 
nels we also take into account the uncertainties in the 
dilution factors due to the imperfect knowledge of the 
levels of the CP-even backgrounds from B~^D^K~, 
D^^K^{K-K-^)non-cj^ and D°^/s:g(7r“7r+7r°)non-A- de¬ 
cays. 

A possible bias in the measured Acp± may come from 
an intrinsic detector charge asymmetry due to asymme¬ 
tries in acceptance or tracking and particle identification 
efficiencies. An upper limit on this bias has been ob¬ 
tained from the measured asymmetries in the processes 
B~^D°h~, and where CP 

violation is expected to be negligible. From the average 
asymmetry, (—1.8±0.9)%, we obtain the limit ±2.7% for 
the bias. This has been added in quadrature to the total 
systematic uncertainty on the CP asymmetry. 

For the branching fraction ratios Rcp± two addi¬ 
tional sources of uncertainty are the correction factors 
used to scale the yield ratios, and the assumption that 
Rcp± = R±/R- The scaling factor, estimated from simu¬ 
lated events, is a double ratio of efficiencies, , 

where denotes the ratio between the selection ef¬ 
ficiencies of and B^D^^p^^Tr. In the 

double ratio the systematic uncertainties arising from 
possible discrepancies between data and simulation are 
negligible, and only the contribution from the limited 
statistics of the simulated samples remains. The as¬ 
sumption Rcp± = R±/R introduces a relative uncer¬ 
tainty ± 2 r^ cos Stt cos 7 on Rcp±, where St^ is the relative 
strong phase between the amplitudes A{B ~and 
A{B ~Since | cos<5^ cosy] < 1 and ± 0.012, 














we assign a relative uncertainty ±2.4% to Rcp±, which 
is completely anti-correlated between Rcp+ and Rcp-■ 
We quote the measurements in terms of x± and 

x+ = —0.082 ± 0.053(stat) ± 0.018(syst), ( 8 ) 

X- = ±0.102 ± 0.062(stat) ± 0.022(syst), (9) 

= -0.12 ± 0.08(stat) ± 0.03(syst). (10) 

The measured values of x± are consistent with those 
found, on a slightly smaller data sample, with the 
B~^DK~ , D—>KgTr~Tr'^ decays, and the precision is 
comparable The measured value of is consistent 
with the upper limits on r from BABAR and Belle Q • 

In conclusion, we have reconstructed B~—s-D^K~ de¬ 
cays with mesons decaying to non-CP, CP-even and 
CP-odd eigenstates. We have improved the previous 
measurements of Rcp± and Acp± [ij, ll^ , and we have 
also expressed the results in terms of the same x± param¬ 
eters as were measured with B~^DK~ , D^KgTr~Tr~^ 
through a Dalitz plot analysis of the D final state 0, 
with a comparable precision. These measurements, com¬ 
bined with the existing measurements of the B^DK de¬ 
cays, will improve the knowledge of the angle 7 and the 
parameter r. 
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